Generation of the second-harmonic Bessel beams via nonlinear Bragg diffraction 
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We generate conical second- harmonic radiation by transverse excitation of a two-dimensional 
annular periodically-poled nonlinear photonic structure with a fundamental Gaussian beam. We 
show that these conical waves are the far-field images of the Bessel beams generated in a crystal by 
parametric frequency conversion assisted by nonlinear Bragg diffraction. 
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The second-harmonic generation, i.e. conversion of two 
photons of the fundamental wave into a single photon 
at twice the frequency, belongs to one of the most in- 
tensively studied nonlinear effects [![. Such a process 
takes place in quadratic optical media, i.e. nonlinear 
media without the center of inversion. It is well estab- 
lished that the efficiency of this as well as other similar 
parametric processes depends critically on the so-called 
phase-matching condition which in birefringent crystals 
is commonly achieved with the temperature or angle tun- 
ing. In media with weak birefringence, an efficient fre- 
quency conversion is achieved by the quasi-phase match- 
ing (QPM) [2[. The QPM technique involves periodic 
modulation of the second-order nonlinearity of the ma- 
terial. In the ferroelectric crystals such as LiNbOs or 
LiTiC>3, this can be easily realized by spatially periodic 
poling [3[. In this way one can create one- or two- 
dimensional nonlinear structures, the so-called pho- 
tonic crystals 0, [Hj]. Depending on the symmetry and 
form of the poling pattern the efficient second-harmonic 
generation (SHG) can be realized for waves propagating 
in particular spatial directions [y]. Typical geometries 
involve, in one spatial dimension, single or multi-period 
rectangular pattern, and, in two spatial dimensions, non- 
linear photonic structures of hexagonal or square symme- 
tries. 

In this Letter, we report on the generation of con- 
ical second-harmonic waves by excitation of a two- 
dimensional annular nonlinear photonic fabricated 
in Stoichiometric Lithium Tantalte (SLT) crystal with a 
single Gaussian beam. We show that the observed rings 
of the second-harmonic radiation actually represent the 
far field of the Bessel beams generated in the crystal via 
the multi-order nonlinear Bragg diffraction. 

We consider a two-dimensional nonlinear photonic 
structure with circular periodic array of ferroelectric 
domains with a constant linear refractive index. In 
general, such two-dimensional QPM photonic structures 
have been explored exclusively in the longitudinal geom- 



etry 0, 0, [I[ of the second-harmonic generation when 
the fundamental beam propagates perpendicular to the 
domains boundaries of periodically varying second-order 
nonlinearity. Here, we explore a novel transverse geom- 
etry when the fundamental light beam propagates along 
the axis of the structure. The front facet of the annular 
poled Stoichiometric Lithium Tantalate sample used in 
our experiments with clearly visible domain boundaries 
is shown in Fig. 1(a) [for technical details, see [7]]. It 
is a Z-cut, L = 0.49 mm thick slab with the QPM pe- 
riod of 7.5 jam and duty factor varying inside the sample 
from 0.7 at Z+ surface to 0.8 at Z-surface. Both, Z+ 
and Z- surfaces are polished. Small 20 nm deep grooves 
that remain after polishing cause weak (less than 3%) 
diffraction of the fundamental beam which we send per- 
pendicular to the domain structure, along its axis. As 
a light source we use a Ti: Sapphire oscillator and regen- 
erative amplifier delivering beam of 140 fs pulses, with a 
repetition rate of 250 kHz and average output power of 
740 mW. Time bandwidth product of the laser pulses is 
0.49. 

The input beam is loosely focused such that the beam 
waist at the input facet of the sample is 147 jam (FWHM) 
which, for a given input power, corresponds to the inten- 
sity about 100GW/cm 2 . The beam covers roughly 10 
domain rings. We observe that the transverse illumina- 
tion of the radial structure by the fundamental beam lead 
to multiple conical emission of the second-harmonic (SH) 
waves [see Figs. 1(b) and 2(a)]. These waves form a set 
of rings in the far field, and their propagation angles sat- 
isfy the Bragg relation, i.e they are determined by the 
ratio of the wavelength of the second harmonic to the 
period of the modulation of the nonlinearity. 

To explain the appearance of the second-harmonic 
rings, we consider the phase-matching conditions Q for 
the corresponding parametric process shown in Fig. [He). 
The general vectorial phase-matching condition can be 
split into two scalar ones: the transverse condition, 
&2 sina m = G m = mGi, and the longitudinal condition, 




FIG. 1: (color online) (a) Front facet of the SLT sample, (b) 
Far-field image of first two diffraction SH rings, (c) Phase 
matching diagram of the SHG process in the transverse QPM 
grating, (d) Polarization structure of the ordinary (left), ex- 
traordinary (middle), and small-angle (right) rings. The fun- 
damental wave (red) has Y polarization. 



&2 cos a m — 2ki = Ak m . With a correct choice of the fun- 
damental wavelength or temperature, a given diffraction 
order will be longitudinally phase-matched (Ak m = 0). 
As an example, in Fig. [Tfc) the third diffraction order 
is phase- matched (Aks = 0), and for this interaction 
both the phase-matching conditions are satisfied simul- 
taneously. Inside the crystal, the propagation angles of 
all transversely phase-matched waves can be found from 
the relation sina m = 171X2 /n2 A, which is analogous to 
the Bragg condition governing the diffraction of waves on 
periodic modulation of the index of refraction. Here the 
refractive index is constant but the nonlinear properties 
of the medium experience spatially periodic variations. 
This effect is known as nonlinear Bragg diffraction, and 
it was first discussed by I. Freund in 1968 [101 ]. 

The condition of the total internal reflection deter- 
mines the maximum number of the second-harmonic 
rings, Mi = A/A2. The SHG process in the transverse 
direction (/3 m = 90°) determines the total maximum 
diffraction orders, M2 = n2^/\2- For this transverse 
SHG process the longitudinal mismatch is Ak m = 2k±, 
and its compensation is possible only when two counter- 
propagating fundamental waves are employed [111 ]. 

For the annular periodically poled structure, the outer 
rings possess interesting polarization properties deter- 
mined by the tensor of the second-order nonlinearity of 
the Lithium Tantalate. This crystal belongs to the 3m 
symmetry point group. For the fundamental beam prop- 
agating along the Z axis, there exist following nonzero 
components: d zxx = d zyy ,d yxx = d xyx = —dy yy . Note 



FIG. 2: (color online) Far field image of experimentally ob- 
served second-harmonic rings. The axis X is horizontal, (a) 
Fundamental beam is exactly in the center of the annular 
structure: first several orders and 13-th order SH ring are 
seen, (b) Dependencies of the wave vector mismatch and the 
predicted SH efficiency normalized to the first-order efficiency 
(in log scale) as a function of the diffraction order. 



that in periodically poled crystals with 3m symmetry the 
values <i xxy , <i yxx , d yyy also change signs periodically [T3 |. 
Two types of the SHG processes are possible: Type I (e- 
00, two ordinary waves generate an extraordinary second- 
harmonic wave) and Type (0-00, two ordinary waves 
generate an ordinary second-harmonic wave). Effective 
nonlinearities for these two processes are [l3[ 



d yyy cos((p + 27), 
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d yyy cos a sm((p + 27) + d zyy sin a, (2) 



where ip is the azimuthal angle measured counterclock- 
wise from the X axis, a is the diffraction angle inside the 
crystal, and 7 is the angle that defines the polarization 
of the input beam (7 = for polarization along X axis) . 
The azimuthal intensity distribution of the "ordinary" 
and "extraordinary" second-harmonic rings is given by 
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where I\ denotes the intensity of the fundamental beam, 
g m = [2/(nm)]sin(7rmD), L cohjm 7r/Ak m and D - rep- 
resents the duty factor. The light intensity in individual 
rings depends also on the phase mismatch Ak m . For 
ring(s) for which both the phase-matching conditions are 
fulfilled the coherence length L co h has to be replaced by 
the actual crystal length L. In Fig. 2(b) we illustrate the 
dependence of the relative intensity (in log scale) in each 
diffraction order on m. 

For small angles a (when tana << 1), the difference in 
the refractive indices of the ordinary and extraordinary 
waves can be neglected, the "ordinary" and "extraordi- 
nary" SH rings perfectly overlap, and as a result the in- 
tensity of conically emitted light does not depend on 7 
and (p. In addition, its polarization becomes linear [Fig. 
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l(d, right)]. For 7 = and for 7 = tt/2 polarization of 
the SH wave coincide with the F-direction (with relevant 
nonlinear components d yxx and d yyy , respectively), while 
for 7 = 7r/4 polarization of the SH wave coincide with the 
X-direction (with relevant nonlinear component d xyx ). 

For relatively large a, the ordinary and extraordi- 
nary waves propagate with different phase velocities 
[n > n e (a)] and at slightly different conical angles. In 
this case, the output polarization depends strongly on 
the actual spatial overlap of the two rings since they 
are orthogonally polarized. In Fig. DJd) we show the 
structure of the predicted polarization on the ring for 
dyyy/d Z yy « —1.7 [14] and a = 20°. For relatively long 
crystals every diffraction order should appear as a dou- 
blet consisting of two orthogonally-polarized rings. 

A full analysis of the SHG process in the QPM struc- 
ture with the azimuthal symmetry shows that, in fact, the 
SH field inside the crystal has a form of a nondiffracting 
Bessel beam [15]. Detailed form depends on the crys- 
tal symmetry. In the considered earlier case of Stron- 
tium Barium Niobate, the SH wave is just the first-order 
Bessel function with a perfect radial polarization [9] . On 
the other hand, in case of LiTa03 or LiNb03, the beam 
structure is more complicated due to the contributions 
of several components of the tensor. However, it 
still can be represented in terms of a superposition of the 
low-order Bessel functions with the angularly modulated 
intensity. The exact formulas are quite cumbersome, but 
they can be simplified in case of small emission angles 
(tana << 1), 

E(p, 0o,z) oc 2nA 2 dyyy Jo(0 [u x sin 27 + u y cos 27] , (4) 

where £ = /c2/?sina and p = (x 2 + y 2 ) 1 ^ 2 is the trans- 
verse radial coordinate. u x and u y are unit vectors along 
X and Y. It is clear that the beam represented by Eq. (|4]) 
is linearly polarized, with its intensity being independent 
of the azimuthal coordinate I oc A-K 2 d 2 yy A A Jq (£). This 
results are in a full agreement with experimental obser- 
vations [Fig. IH[d)]. 

The observed SH rings for the fundamental wavelength 
of 822 nm are shown in Fig. EJa). The measured conical 
angles agree well with those determined from the non- 
linear Bragg diffraction formula. By re-scaling the size 
of the lower-order diffraction pattern we found that the 
phase matched SH emission which is represented in the 
far-zone by the ring with the largest diameter and coni- 
cal angle of 45.6° corresponds to the high 13-th diffrac- 
tion order. In Fig. EJb) we depict the calculated (us- 
ing the published refractive index data [16]) longitudinal 
phase mismatch and relative power of the SH waves for 
A = 836 nm. The highest peak corresponds to the phase- 
matched (Afci3 = 0) 13-th order SH ring. 

The phase-matching curve is shown in Fig. [3^a) to- 
gether with the theoretical prediction based on the model 
from Ref. [§]. We observe good qualitative agreement 
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FIG. 3: (color online) (a) Experimentally measured (circles) 
and theoretically predicted (solid line) phase matching curves, 
(normalized to the maximum signal of the experiment), (b) 
Quadratic dependence of the 13-th order SH ring on input 
power, (c) Spectrum of the fundamental pulse, (d) Spectra 
of the SH rings for non phase matched (first-order) and phase 
matched (13-th order) SH ring. 



between the experimental curve and the theoretical pre- 
diction. In Fig. [3](b) we verify a quadratic dependence of 
the SH power of the phase matched 13-th SH ring as a 
function of the input power. From this figure (obtained 
at 640 mW fundamental average power taking into ac- 
count all losses), we obtain the internal efficiency for the 
phase-matched emission of 7713 = 0.0056% / W. Such low 
efficiency is attributed to the fact that we use high order 
QPM interaction. Additionally, the variations of the duty 
factor inside the crystal leads to an effective length 2-3 
times less than the real sample thickness. The measured 
efficiency of the SH power generated in the first-order 
ring is 220 times less than 7713. The spectra are shown 
in Fig. 3(c,d). The fundamental spectral width is 8nm. 
The spectra of the SH ring in the case of phase-matching 
of the 13-th order is 1.3 nm; i.e. almost 2 times less than 
expected for a stationary process. This indicates that 
the SHG process takes place in the group-velocity mis- 
match conditions. Indeed, the calculated group-velocity- 
mismatch length is 0.092 mm for 140 fs fundamental pulse 
width, that is a half of the expected effective length of the 
sample. Subsequently, the output spectra (pulse length) 
is expected to be reduced (stretched) by a factor of 2. In 
contrast, the spectra of non-phase matched rings [see Fig. 
3(d)] do not show any reduction and are close to those 
predicted for the stationary process, Ac^sh = V^Ac^ff- 

The rings are sensitive to the axial alignment of the 
fundamental beam. With the off-center beam, the rings 
transform into a two-peak structure, the efficiency in- 
creases, and the number of the observed diffracted or- 
ders increases too. For a horizontal shift of the sample, 
the diffraction spots appear horizontally. For a vertical 
shift of the sample, the diffraction spots appear vertically. 
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FIG. 4: (color online) (a-c) Modification of the diffraction 
pattern as a function of the horizontal misalignment of the in- 
cident beam. Numbers indicate the shift in /mi. The cry st al- 
lographs X axis is horizontal. (d,e) Experimentally recorded 
intensity distribution of m — 1 and m — 13 SH rings, (f) The- 
oretical prediction for the azimuthal intensity distribution of 
the m — 13 SH ring evaluated from Eq. (3). 




Analyzer angle 



Input polariser angle 



FIG. 5: (color online) Polarization properties of the SH rings: 
(a) Power of the 1-st order SH ring vs. analyzer angle, (b) 
Total power of the 1-st and 13-th order SH rings vs. the input 
polarization angle. 



This is explained by the fact that the structure resem- 
bles a one-dimensional grating formed from the sectors 
situated close to the horizontal (vertical) diameter of the 
rings. The modifications of the diffraction pattern for the 
phase-matched 13-th order ring in the case of off-center 
excitement are shown in Figs.[4](a-c). Figures [H(d,e) show 
the experimental intensity distribution of the 1-st and 13- 
th order rings. The theoretical prediction for the 13-th 
order ring is given by Eq. (J3]), and it is shown in Fig.[U[f). 

As follows from Fig. 1(b) and Figs. [U(a,d,e), the 
second-harmonic rings exhibit an additional azimuthal 
modulation with 6 well defined peaks. These peaks are 
insensitive to the input polarization and intensity, and 
they appear even in the regime of linear diffraction of 
the fundamental beam on the 20 nm deep surface relief - 
replica of poling pattern. The peaks are artifacts of the 
discrete character of domain boundaries resulting from 
the well known tendency of ferroelectric domains in the 



SLT to retain their hexagonal shape @, QjJ . 

In Fig. 5 we show our results of the polarization mea- 
surements of the second-harmonic rings. For the linearly 
polarized input fundamental beam the generated wave 
is also linearly polarized in the case of the first-order 
ring. This result agrees with Eqs. (1) and (2) in the 
limit a — > 0. Polarization state becomes more complex 
for higher-order nonlinear diffraction as the output signal 
contains both ordinary and extraordinary contributions. 

In conclusion, we have reported on the first observation 
of the second-harmonic Bessel beams generated by two- 
dimensional annular nonlinear photonic structures. We 
have explained the effects observed in experiment by em- 
ploying the concept of nonlinear Bragg diffraction com- 
bined with the longitudinal and transverse phase match- 
ing conditions. We have studied the polarization prop- 
erties of the conical second-harmonic radiation and de- 
scribe how the beam polarization depends on the crystal 
symmetry and conical angle. 
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